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Abstract 

The present work reports results for: i) pd radiative capture observables 
measured at center-of-mass (cm.) energies in the range 0-100 keV and at 
2 MeV by the TUNL and Wisconsin groups, respectively; ii) contributions 
to the Gerasimov-Drell-Hearn (GDH) integral in ^He from the two- up to 

1 



the three-body breakup thresholds, compared to experimental determinations 
by the TUNL group in this threshold region; iii) longitudinal, transverse, 
and interference response functions measured in inclusive polarized electron 
scattering off polarized ^He at excitation energies below the threshold for 
breakup into ppn, compared to unpolarized longitudinal and transverse data 
from the Saskatoon group. The calculations are based on pair-correlated- 
hyperspherical-harmonics bound and continuum wave functions obtained from 
a realistic Hamiltonian consisting of the Argonne vis two-nucleon and Urbana 
IX three-nucleon interactions. The electromagnetic current operator includes 
one- and two-body components, leading terms of which are constructed from 
the Argonne vig interaction (specifically, its charge-independent part). Two- 
body currents associated with A-isobar degrees of freedom are treated non- 
perturbatively via the transition-correlation-operator method. The theoret- 
ical predictions obtained by including only one-body currents are in violent 
disagreement with data. These differences between theory and experiment 
are, to a large extent, removed when two-body currents are taken into ac- 
count, although some rather large discrepancies remain in the cm. energy 
range 0-100 keV, particularly for the pd differential cross section a{6) and 
tensor analyzing power T2q{9) at small angles, and contributions to the GDH 
integral. A rather detailed analysis indicates that these discrepancies have, 
in large part, a common origin, and can be traced back to an excess strength 
obtained in the theoretical calculation of the Ei reduced matrix element as- 
sociated with the pd channel having L, S, J = 1,1/2,3/2. It is suggested that 
this lack of Ei strength observed experimentally might have implications for 
the nuclear interaction at very low energies. Finally, the validity of the long- 
wavelength approximation for electric dipole transitions is discussed. 
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I. INTRODUCTION 



Radiative capture, photo- and electro-disintegration reactions are a useful tool for explor- 
ing the structure of nuclei and their electromagnetic response. The theoretical description 
of these processes requires knowledge of the nuclear bound- and scattering-state wave func- 
tions and electromagnetic transition operators. In this respect, the trinucleon systems play a 
unique role because of the capability, achieved in the last few years, to obtain very accurate 
wave functions for both bound and continuum states from realistic Hamiltonian models . 

The accuracy of the calculated trinucleon continuum wave functions has been verified 
by comparing results for a variety of Nd scattering observables obtained by a number of 
groups using different techniques (see Ref. ^ and references therein). In fact, good overall 
agreement exists between theory and experiment for both elastic and inelastic Nd cross 
sections and polarization observables, with the only notable exceptions of the pd and nd 
vector analyzing powers at low energies, which are both underpredicted by theory at the 
30% level Indeed, the Ay "puzzle" constitutes an excellent example of how, once 

the numerical uncertainties in the calculation of the continuum wave functions have been 
drastically reduced, Nd scattering observables can be used to study the sentitivity to two- 
and three-nucleon interaction models. 

Electromagnetic processes, such as the low-energy pd radiative fusion and threshold 
photo- and electro-disintegration of '^He under consideration in the present work, provide 
additional and essential insights into the structure of the trinucleons, since they can be 
used to test and refine (or possibly, discriminate among) models for the nuclear interactions 
and electromagnetic transition operators. Furthermore, they allow us to study a number of 
other, closely related issues. A specific example of these is the relevance of A-isobar degrees 
of freedom for the proper description of photo- and electro- nuclear observables [^]. 

There are now available many high-quality data, including differential cross sections, 
vector and tensor analyzing powers, and photon polarization coefficients, on the pd radiative 
capture at cm. energies ranging from to 2 MeV The goal of the present study is 
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to determine the extent to which this rich body of data can be described satisfactorily by a 
calculation based on a realistic Hamiltonian (consisting of the Argonne fig two-nucleon ||T^ 
and Urbana-IX three- nucleon |Tl| interactions), and a current operator including one- and 
two-body components, leading terms of which are constructed consistently with the two- 
nucleon interaction. The present work updates and extends previous ones, which only dealt 
with pd radiative capture in the cm. energy range 0-100 keV and nd radiative capture at 
thermal neutron energies It contains: i) an improved treatment of the pd continuum 
within the pair-correlated-hyperspherical-harmonics (PHH) scheme; ii) calculations of pd 
radiative capture polarization observables in the cm. energy range 0-100 keV, measured by 
the TUNL group 0, and at 2 MeV, measured by the Wisconsin group 0; iii) a calculation 
of the contribution to the Gerasimov-Drell-Hearn integral of ^He in the threshold region. 



compared to results of an experimental determination by the TUNL group [|12|; and iv) a 
calculation of the longitudinal, transverse, and interference response functions measured in 
inclusive scattering of polarized electrons off polarized ^He for excitation energies below the 
three-body breakup threshold and for momentum transfers up to 5 fm~^, compared with 



longitudinal and transverse data measured by Retzlaff et al. ||T^ at g=0.88, 1.54 and 2.47 
fm~^. 

The paper is organized into five sections and an appendix. In Sec. II we discuss the 
calculation of the bound and scattering wave functions with the PHH method, and summa- 
rize a number of results obtained for Nd elastic scattering observables, comparing them to 
experimental data. In Sec. Ill we briefly review the model for the electromagnetic transition 
operator, while in Sec. IV we present an exhaustive comparison between theory and exper- 
iment for all available pd capture, photo- and electro-disintegration data in the threshold 
region. Finally, in Sec. V we summarize our conclusions. A collection of formulas for the cal- 
culation of various observables from the reduced matrix elements of the current and charge 
operators, which complement those published earlier in Ref. [0], are given in the Appendix. 
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II. BOUND- AND SCATTERING-STATE WAVE FUNCTIONS 



The ^He bound-state and pd scattering-state wave functions are obtained variationally 
with the pair-correlated-hyperspherical-harmonics (PHH) method from a reahstic Hamil- 
tonian including the Argonne vis two-nucleon [10|| and Urbana-IX three-nucleon |11 



in- 



teractions (the AV18/UIX model). The PHH method, as implemented in the calculations 
reported in the present work, has been developed by Kievsky, Rosati and Viviani in a series 
of papers appeared in various journals between 1993 and 1995 PJ^JT^. Here it will be re- 
viewed briefly for completeness, and a summary of relevant results obtained with it for the 
three-nucleon bound-state properties and Nd scattering observables at energies below the 
three-body breakup thresholds will be presented. 



The three-nucleon bound-state wave function \l/3 is expanded as [14 



^3 = E^^ E ZM^;Jk) , (2.1) 

a,K P cyclic ijk 



where p is the hyperradius, p = yx^+yf with the Jacobi variables Xj and defined, 
respectively, as Xj = rj — and = (r^ + — 2vi)/^/3. The known functions ZaK{i]jk) 
are antisymmetric under the exchange j ^ k and account for the angle-spin-isospin and 
hyperangle dependence of channel a, K. The index a denotes collectively the spectator i 
and pair jk orbital and spin angular momenta and isospins coupled to produce a state with 
total angular momentum and isospin J'", T = |, while the index K specifies the order of 
the Jacobi polynomial in the hyperangle cos0j = Xi/p. Correlation factors, which account 
for the strong state-dependent correlations induced by the nucleon-nucleon interaction, are 
included in the functions ZaKii'^jk). 

The Rayleigh-Ritz variational principle, 

(5„*3|i?-^3|^3) = , (2.2) 

is used to determine the functions Uaxip)- Carrying out the variations with respect to the 
UaK^ leads to a set of coupled second-order differential equations. After discretization in 
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the variable p, this set of differential equations is converted into a generalized eigenvalue 
problem, which is then solved by standard numerical techniques P |. 



Fully converged AV18/U1X PHH wave functions predict three-nucleon binding energies 
and matter radii respectively given by Bi^R) =8.48 MeV and {r^i^R)^/^ = 1.725 fm, 
^(^He) = 7.73 MeV and {r'^i^Re)^/'^ = 1.928 fm, in agreement with corresponding "exact" 
Green's function Monte Carlo results and with the available experimental data. 

The Nd cluster wave function "^1^2"^% having incoming orbital angular momentum L 
and channel spin S* (S* = 1/2, 3/2) coupled to total JJz, is expressed as 

^LSJJ. ^ ^JJ. ^ ^LSJJ, ^ (^2.3) 

where vanishes in the limit of large intercluster separation, and hence describes the 
system in the region where the particles are close to each other and their mutual interactions 
are large. In the asymptotic region, where intercluster interactions are neglig ible, ^^^^^^ 
(for pd, as an example) is written as 



cyclic ijk L' S' 



X 



Sll'Sss' + RLS,L's'iP) -^girpd) 

prpd prpd 



(2.4) 



where (pd is the deuteron wave function, p the magnitude of the relative momentum between 
deuteron and proton, and Fl and Gl are the regular and irregular Coulomb functions, 
respectively. Note that for nd scattering, Fi{x)/x and Gl{x)/x are to be replaced by the 
regular and irregular spherical Bessel functions. The function gij-pd) modifies the Giiprpd) 
at small Vpd by regularizing it at the origin, and girpd) — 1 as Vpd > 10 fm, thus not 
affecting the asymptotic behavior of \E'f+2"^^- Finally, the real parameters R'ls us'ip) 
i?-matrix elements which determine phase-shifts and (for coupled channels) mixing angles 
at the energy j9^/(2/i), /x being the 1+2 reduced mass. Of course, the sum over L'S' is over 
all values compatible with a given J and parity. 

The "core" wave function \E'c' is expanded in the same PHH basis as the bound-state 
wave function \E'3, Eq. ( |2.1[ ), and both the matrix elements Rls L'S'iP) ^^"^ functions Uaxip) 



occurring in the expansion of are determined by making the functional p|Q 

[Ris,L'S'ip)] = Ris,L'S'ip) - (^S'-^'li^ -E,- l^l^f+T^) , (2.5) 

stationary with respect to variations in the Rlsl's' ^^"^ '^aK (Kohn variational principle). 
Here E2 = —2.225 MeV is the deuteron energy. 

Phase-shifts and mixing angles for nd scattering at energies below the three-body breakup 
threshold have been obtained from a realistic Hamiltonian model, and have been shown 



to be in excellent agreement with corresponding Faddeev results |]15[, thus establishing 
the high accuracy of the PHH expansion for the scattering problem. It is important to 
emphasize that the PHH scheme, in contrast to momentum-space Faddeev methods, permits 
the straightforward inclusion of Coulomb distortion effects in the pd channel. The PHH 
results for pd elastic scattering are as accurate as those for nd scattering. 

The nd and pd doublet and quartet scattering lengths predicted by the AV18 /UIX model 
are listed in Table |, and are found to be in excellent agreement with the available experi- 
mental values. 



III. THE ELECTROMAGNETIC CURRENT OPERATOR 

The nuclear charge and current operators consist of one- and two-body terms that operate 
on the nucleon degrees of freedom: 

p(q) = Epl'^(q) + Epl?(q) > (s-i) 

i i<j 

j(q) = Ejr^(q) + EjS?(q) > (3-2) 

i i<j 

where q is the momentum transfer. 

The one-body operators p\^^ and j-^"* have the standard expressions obtained from a 
relativistic reduction of the covariant single-nucleon current, and are listed below for conve- 
nience. The charge operator is written as 

pf^(q) = pSR(q) + pSc(q) , (3-3) 
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with 



ei e 



iqri 



(3.4) 



p!,'iS,c(q) 



, , - 1 e^'^-' - (2 p, - 6,) q ■ (CT, X p,) e'^^" , (3.5) 



where g„ = q —u is the four-momentum transfer, and u is the energy transfer. The current 



operator is expressed as 



1 



(3.6) 



where [■ ■ ■ , ■ ■ ■]+ denotes the anticommutator. The following definitions have been intro- 
duced: 



1 

1 



GlMD + GlMDr. 



(3.7) 
(3.8) 



and p, cr, and r are the nucleon's momentum, Pauli spin and isospin operators, respectively. 
The two terms proportional to 1/ m? in Pj- ^^^^ known Darwin- Foldy and spin-orbit 
relativistic corrections |T^,|T^, respectively. 

The dipole parametrization is used for the isoscalar {S) and isovector iV) combina- 
tions of the electric and magnetic nucleon form factors (including the Galster form for the 



electric neutron form factor [T^). It is worth emphasizing that the neutron form factors, 
particularly the electric one, are not well known and, therefore, the available semi-empirical 
parameterizations for them differ widely, particularly at high momentum transfers. Until 
this uncertainty in the detailed behavior of the electromagnetic form factors of the nucleon is 
narrowed, quantitative predictions of electronuclear observables at high momentum transfers 
will remain somewhat tentative. 

The most important features of the two-body parts of the electromagnetic current oper- 
ator are summarized below. The reader is referred to Refs. P,p^ for a derivation and listing 
of their explicit expressions. 
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A. Two-body current operators 



The two-body current operator has "model-independent" and "model-dependent "compo- 
nents, in the classification scheme of Riska [^]. The model- independent terms are obtained 
from the charge-independent part of the AV18, and by construction satisfy current 
conservation with this interaction. The leading operator is the isovector "vr-like "current ob- 
tained from the isospin-dependent spin-spin and tensor interactions. The latter also gener- 
ate an isovector "p-like "current, while additional model-independent isoscalar and isovector 
currents arise from the isospin-independent and isospin-dependent central and momentum- 
dependent interactions. These currents are short-ranged and numerically far less important 
than the vr-like current. 

The model-dependent currents are purely transverse and therefore cannot be directly 
linked to the underlying two-nucleon interaction. The present calculation includes the 
isoscalar p7r7 and isovector co'7r7 transition currents as well as the isovector current asso- 
ciated with excitation of intermediate A-isobar resonances. The piT'-f and uott'-) couplings are 



known from the measured widths of the radiative decays p ^ ti'J p2| and — » 7r7 p3 
respectively, while their momentum-transfer dependence is modeled using vector-meson- 
dominance. Monopole form factors are introduced at the meson-baryon vertices with cutoff 
values of A^=3.8 fm^^ and Ap=A^=6.3 fm^^ at the ttNN, pNN and ujNN vertices, respec- 
tively. 

Among the model-dependent currents, those associated with the A-isobar are the most 
important ones. In the present calculation, these currents are treated within the transition- 



correlation-operator (TCO) scheme developed in Ref. [^. In such an approach, the A 



degrees of freedom are explicitly included in the nuclear wave functions by writing 



Af+A — 



jTR 
ij 

i<j 



^ , (3.9) 



where \& is the purely nucleonic component, S is the symmetrizer and the transition cor- 
relations U^j^ are short-range operators, that convert NN pairs into A^A and AA pairs. 
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In the present study the \I/ is taken from PHH solutions of the AV18/UIX Hamiltonian 
with nucleons only interactions, while the Ufj^ is obtained from two-body bound and low- 
energy scattering state solutions of the full A^-A coupled-channel problem. This aspect of 
the present calculations as well as the justification for going beyond the traditional pertur- 
bative treatment of A degrees of freedom in nuclei, were discussed at length in the original 



work p4[, and have been reviewed most recently in Ref. p5[, making a further review here 



unnecessary. 



In the TCO approach, both 7A^A and 7AA Mi couplings are considered |£4|. The 



values used for these couplings are, respectively, /x^ata = 3 n.m. and /x^-aa = 4.35 n.m. The 



former is taken from an analysis of 7A^ data in the A-resonance region p6[ , while the latter 
is obtained from a soft-photon analysis of pion-proton bremsstrahlung data near the A^^ 



resonance 



Electromagnetic observables require evaluation of matrix elements of the type [p^,p5 





j 








^A^+A.)]^/^ 



(3.10) 

where the wave functions and currents include both nucleonic and A-isobar degrees of free- 
dom. To evaluate such a matrix element, it is convenient to expand the wave function ^I/at+a 

as 



^7V+A = ^ + Ef^r^ + -- - ' (3-11) 

i<j 

and write the numerator of Eq. (|3.10|) , in a schematic notation, as 



(^Tv+Aj I j I ^^+A,i) = (^/ I j(Aronly) | + | j(A) | ^<>) , (3.12) 

where j(A^only) denotes all one- and two-body contributions to j(q) which only involve 
nucleon degrees of freedom, i.e., j(A^only) = j(^)(A^ ^ A^) + p\NN A^A^). The 
operator j(A) includes terms involving the A-isobar degrees of freedom, associated with the 
explicit A currents j''"'^^(A^ A), j^^''(A N), and j^^''(A A), and with the transition 
operators Ulj^. Of course, the presence of A admixtures in the wave functions also influences 
their normalization. 
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Finally, we note that the contributions associated with the p7r7, Ci;7r7 and A-excitation 
mechanisms are, at low and moderate values of momentum transfer (g < 5 fm~^), typically 
much smaller than those due to the leading model-independent 7r-like current P,|25|. 

B. Two-body charge operators 

While the main parts of the two-body currents are linked to the form of the two-nucleon 
interaction through the continuity equation, the most important two-body charge operators 
are model-dependent, and should be considered as relativistic corrections. Indeed, a consis- 
tent calculation of two-body charge effects in nuclei would require the inclusion of relativistic 
effects in both the interaction models and nuclear wave functions. Such a program is just at 
its inception for systems with A> 3. There are nevertheless rather clear indications for the 
relevance of two-body charge operators from the failure of the impulse approximation (lA) 
in predicting the deuteron tensor polarization observable ||28[, and charge form factors of the 
three- and four-nucleon systems p9| , |25[| . The model commonly used includes the vr-, p-, 
and tu-meson exchange charge operators with both isoscalar and isovector components, as 
well as the (isoscalar) pn'y and (isovector) uTC'-f charge transition couplings, in addition to 
the single-nucleon Darwin-Foldy and spin-orbit relativistic corrections. The vr- and p-meson 
exchange charge operators are constructed from the AV18 isospin-dependent spin-spin and 
tensor interactions, using the same prescription adopted for the corresponding current op- 
erators. It should be emphasized, however, that for g<5 fm~^ the contribution due to the 
TT-exchange charge operator is typically an order of magnitude larger than that of any of the 
remaining two-body mechanisms and one-body relativistic corrections. 

IV. RESULTS 

The present section contains results for the pd radiative capture in the cm. energy range 
0-100 keV, compared to the TUNL data and at 2 MeV, compared to the Wisconsin 
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data , as well as results for the threshold electrodisintegration of ^He at momentum trans- 
fers g=0.88, 1.64 and 2.74 fm~^, compared to Bates data [0. It also contains predictions for 
the contribution to the Gerasimov-Drell-Hearn integral of '^He up to the three-body breakup 
threshold, compared to a very recent experimental determination of this contribution in the 
pd threshold region by the TUNL group . 

The calculations are based on bound- and scattering-state PHH wave functions, obtained 
from the Argonne f is two-nucleon and Urbana-IX three-nucleon interactions (see discussion 
in Sec. H). The model for the electromagnetic current operator includes one- and two-body 
components, leading terms of which are constructed from the two-nucleon interaction (in 
the present case, the charge independent part of the AV18). The currents associated with 
the excitation of A-isobars are treated non-perturbatively within the transition-correlation- 



operator scheme, as sketched in Sec. |Tl. 

A study of the pd capture cross section and polarization observables in the cm. energy 
range 0-100 keV was reported earlier in Ref. 0. In the present work, improvements in the 
PHH variational treatment of the P-wave pd channel have led to significant changes in the 
predictions for the S'-factor and some of the polarization observables, particularly the vector 
analyzing power Ay, than previously published 0. We therefore provide an update of that 
study. 



A. The pd capture below 100 keV cm. energy 

The observed linear dependence upon the energy of the S'-factor and the observed angular 
distributions of the polarization observables indicate that the pd radiative fusion proceeds, 
at these low energies, through S- and P-wave captures. Therefore, the contributing reduced 

oil 0-- 0-- 1-- 1-- 

matrix elements (RMEs) are M-^ ^ ^ , ^ ^ , and -E2 ^ ^ S-wave capture, and Ei^^ , E^^'^ , 

ill ^33 -|^1.3 -j^33 

E^^^ , Ei^^ , M2 ^ ^ , and ^ ^ in P-wave capture, where and Ei are the magnetic and 
electric £th-pole operators, respectively. The superscripts LSJ refer to the relative orbital 
angular momentum L between the p and d clusters, the channel spin S {S = 1/2 or 3/2), 
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and the total angular momentum J (J = L + S), respectively. The M2 transition from the 
L, 5, J =1,3/2,5/2 capture state has been neglected. 

Improvements in the PHH description of the pd continuum wave functions, particularly 
in the short-range part, have led to significative changes in the values of the P-wave cap- 
ture RMEs than previously calculated [Q. The S-wave capture RME values have remained 
essentially unchanged. 

The El RMEs have in fact been calculated in two different ways: firstly, by direct 
evaluation of the matrix elements of the Ei multipole operator, 

^^^^^-(CTI^^^3^^^^^ (4.1) 

Ei, = - [ f/x j(x) ■ V X ji(gx)Y|i(x) , (4.2) 
q J 

where j(x) is the nuclear current density operator, ji{qx) is the spherical Bessel function 
of order one, Y{j^(x) are vector spherical harmonic functions, and the pd wave function is 
constructed to satisfy outgoing- wave boundary conditions as in Eq. (4.3) of Ref. 0. 

Secondly, in the long- wavelength- approximation (LWA) (certainly justified in the energy 
range under consideration here) the Ei operator can be expressed as 



Eia^^ia(LWA1) + Eia(LWA2) + Eia(LWA3) , (4.3) 



where 



/2 

Eu(LWAl) = 



H , J dxxyiA(x) p(x) 



(4.4) 



2 

EiA(LWA2) = ^y'rfxxFiA(i)x-j(x) , (4.5) 



3^2 

Eia(LWA3) 



H, Jd^x^Yix{k)p{^) 



(4.6) 



15 

Here the continuity equation has been used to relate V ■ j(x) occurring in _Eia(LWA1) and 
i?iA(LWA3) to the commutator —i[H, p(x)], where p(x) is the charge density operator. 
Evaluating the RMEs of these operators leads to 
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^LSJ _ E^^-^{LWAl) + ^*^(LWA2) , (4.7) 

with 

Ef^^(LWAl) = JJj^^§^^(n'''\ J rfxxn,(x)p(x)|^f^f ^) , (4.8) 



Et''{LWA2) = / d^xY,,i±)^-j{^)\^'4'/^) , (4.9) 

where only terms up to order have been retained, and it has been assumed that the 
initial and final PHH wave functions are exact eigenfunctions of the Hamiltonian, so that 
[if, p(x)] —qp{x) in the matrix element. The contribution to the Ei RMEs associated 
with the LWA3 operator, defined in Eq. ([4.6|) , is of order q^, proviso the assumption above. 

By ignoring two-body contributions to the charge density operator, we further approxi- 
mate p(x) ~ pjjp^(x) + pj^Q(x), and write correspondingly 

^^'■^■^(LWAl) ~ Et^-^ilMAc) + Ei^'^-^(LWAb) , (4.10) 

where 

Ef-^(LWAc) = f , 4^(^l"lE^.<Al^f'2'^) , (4.11) 

Hence, up to order g^, the Ei RMEs in LWA are given by Eq. (|47|) , with £'f'^'^(LWAl) 
defined in Eqs. ( [4.10| )- (|4.12| ). We re-emphasize that the currents used in the present work 



satisfy, by construction, the continuity equation with the AV18 interaction, namely q-j(q) = 
T + vis, PNR(q)l • Therefore, if the contributions Ef'^-^(LWAb) and ■^■^(LWA2) (the latter 
of order g^) were to be negligible, the degree of agreement between the £'i(LWAc) and 
full current results, obtained directly from Eq. ( [4.1| ), would simply refiect the extent to 
which the present variational wave functions are truly exact eigenfunctions of the AV18/UIX 
Hamiltonian. In Table || the results for the Ei RMEs obtained by direct evaluation of the 
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El multipole operator matrix elements given in Eq. are compared with those obtained 
in LWA. Note that the ii^;f"^'^(LWA2) contribution, Eq. ( [4.5| ), has been estimated here by 
using only the spin-part of the lA current operator. Of course, single-nucleon convection as 
well as two-body currents provide additional corrections to j(x), but these have been ignored 
in the evaluation of Eq. ( |4.9| ). Finally, note that the RMEs listed in Table [I^ are related to 
those defined in Eqs. Q and (PD-dD via 



where /iat is the nuclear magneton and Vj-ei is the pd relative velocity. The quantities X^J"^ 
are easily shown to remain finite in the limit frei —>■ 0. 

As can be seen by inspection of Table for the S = 1/2 states the values of the LWAc 
and "full" El RMEs are very close to each other and the remaining small differences between 
them are presumably due to the approximate calculation of the E'f '^■^(LWA2) contribution 
carried out here (see above). Therefore, the present PHH wave functions for these P-wave 
channels are a good approximation to the true eigenfunctions. This is to be contrasted 
with the results reported in Table VI of Ref. [0, where significant differences remained 
particularly for the S =1/2 J = 1/2 Ei RME. 

The situation is different for the S = 3/2 states. Here, as discussed in Ref. ^ and below, 
due to cancellations among different contributions, the ii^i(LWAc) RMEs are rather small 
and indeed have similar magnitude as the i?i(LWAb) and -Ei(LWA2) RMEs. Hence, if the 
PHH wave functions for these channels also approximate well the true eigenfunctions-and 
indeed there is no reason to believe that this is not so-, then the differences between the 
"full" and LWAc values provide a "measure" of the corrections beyond the standard LWAc. 
Inspection of Table |I| in fact shows that the leading LWAc form of the Ei operator is 
inadequate for these channels. This point will be further elaborated below and in the next 
subsection. 

As a final remark, note that the continuity equation requires the presence of three-body 
currents associated with the three-nucleon interaction. These currents have been studied in 




(4.13) 
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Ref. p5|], where they were found to give a very small contribution to the trinucleon magnetic 
moments and form factors. It is, therefore, unlikely that three-body current contributions 
influence significantly the "full" predictions discussed above, as it was speculated in Ref. 0. 

The leading RMEs are the doublet and quartet Mi in S-wave capture and doublet 
(namely, S = 1/2) Eis with J = 1/2 and 3/2 in P-wave capture. The Mi and doublet 
El strengths are comparable. The E2 RME is more than an order of magnitude smaller 
than any of the two Mi's, the quartet (namely, S = 3/2) Eis with J = 1/2 and 3/2 are 
an order of magnitude smaller than the doublet Eis, and the M2 strength is negligible. In 
LWAc the ii^i-multipole operator is spin- independent, and transitions from the ^He ground 
state to the S=3/2 channel pd-states are inhibited, since they must proceed through the 
relatively small D-wave components of the ^He wave function. Hence, the quartet Ei RMEs 
are individually small. Such is not the case for the doublet Ei RMEs, which result from 
transitions involving the S'=l/2 pd states and the dominant S-wave component of the ^He 
ground state. 

The effects due to two-body currents and A-degrees of freedom are large on the doublet 
Ml RME-at Ep = 40 keV, they increase it in magnitude by 87%-and significant on the 
quartet Mi RME^at Ep = 40 keV, they reduce it by 10%. The eI^^ and e}^^ RMEs are 
increased by 7% and 10%, respectively, by two-body current contributions (at Ep = 40 keV); 
however, these contributions dominate the quartet Ei^ , interfering destructively with the 
one-body (lA) results, see Table 0. 

The calculated S'-factor is compared with the TUNL data J^,^ in the lab energy range 
Ep=0-300 keV in Fig. |l[ The present results are in better agreement with data for Ep 
between 40 and 80 keV than previously reported 0. However, at higher energies the calcu- 
lation is slightly above the data, in particular for the data point at Ep=160 keV. Below 40 
keV theory is again slightly above the data of Ref. [0]. 

The calculated S'-factor value in S-wave capture at zero energy is 0.110 eV-b, in excellent 
agreement with the corresponding experimental value 0.109 ± 0.010 eV-b However, 
the calculated value in P-wave capture is 0.109 eV-b, which substantially overpredicts the 
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experimental value 0.074 ± 0.01 eV-b 0. This overprediction at zero energy (and below 40 
keV) is due to excess strength in the calculated L, S, J = 1, 1/2, 3/2 Ei RME, as discussed 
below. 

The predicted angular distributions of the differential cross section cr{6), vector and ten- 
sor analyzing powers Ay{6) and T2q{6), and photon linear polarization coefficient Pry{6) are 
compared with the TUNL data from Ref . in Fig As in Ref . , we have integrated 
the theoretical results, weighted with the energy dependence of the cross section and target 



thickness, for the purpose of comparing them with experiment ||3TI|. The ~ 10% changes 
in the Ei RMEs, due to the use of more accurate PHH wave functions, are responsible for 
the improved description of the vector analyzing power Ay which was significantly under- 
predicted by theory in our earlier work 0. However, the calculated a (9) {T2o{9)) is much 
smaller (larger) than the experimental values at small angles. The small angle discrepancy 
for T20 was also present in Ref. 0, although it was not as pronounced as found here. 

Recently, the TUNL group has measured additional polarization observables at Ep=AO 



keV [0, specifically the tensor analyzing powers iTn, T21, T22 and the circular polarization 
asymmetry coefficient A^. All these observables are found to be in satisfactory agreement 
with the present calculations [O. 

The unpolarized cross section a (6) and the tensor analyzing power are given by 

a{e) = Y.a,P,{cose) , a{e)T2^{e) = J2cePe{cose) , (4.14) 
e=o £=0 

where the Pe are Legendre polynomials, and the coefficients ae and q can be expressed in 

terms of the various RMEs 0. Hereafter, for ease of presentation, we introduce the notation: 

m2j+i = M','' , (4.15) 
P2J+1 = eI^' , (4.16) 
q2j+i = e\^' , (4.17) 

the energy dependence of these RMEs being understood. The leading coefficients in the 
expansion of T20 are cq and C2, and their expressions are: 



Co = CTl 



2^ip2q;)-J-\q2\'-J-'fi{p4ql) 



ml 



(4.18) 



C2 = (Ti 



-3?(m2m4) + 
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(4.19) 



where 
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(4.20) 



Note that the expression for cq reported in Eq. (B7) of Ref. contains several misprints 
(the different expression for cti used in Ref. accounts for the fact that there the Legendre 
coefficients were given in terms of the RMEs X/"^'^ defined in Eq. (6.3) of that work). 

The coefficients cq and C2 are found to be rather sensitive to the q2 and q^ RMEs (in par- 
ticular, Co vanishes if both q2 and q^ are set to zero). As discussed above, these RMEs result 
from transitions connecting small components of the wave functions, and are dominated by 
many-body current contributions (see Table |I|). It is not clear, at this stage, whether the 
discrepancies between the measured and calculated T20 (particularly pronounced at small 
angles) are due to deficiencies in the wave functions, or rather the interactions generating 
the wave functions, and/or the many-body current models. 

The expressions for the leading coefficients Oo and 02 in the Legendre expansion of a (6) 

are 



ao = ai [|m2p + |m4p + |p2p + |g2p + 1^4^ + |g4p 
a2 = (Ti 



V2^iP2Pl) - 1\P,\' - -^^{q2Pl) + l\ 



(4.21) 
(4.22) 



For these coefficients the contributions of the q2 and ^4 RMEs are completely negligible, 
and the observed small-angle discrepancy between theory and experiment is due to excess 
strength in the calculated p^ RME, see below. In particular, we have verified that in this 
low-energy regime the differences between the "full" results and the results obtained by using 



19 



the LWAc values for the P2, Pa, (I2 and RMEs are small for all the observables considered in 
Fig. 1^. This is not the case for the pd capture at 2 MeV, as discussed in the next subsection. 

The extensive body of data measured at TUNL has allowed the determination of the lead- 
ing Ml and Ei RMEs (magnitudes and phases) via fits to the measured observables p!2| , |3T| . 
The results of this fitting procedure are compared with the calculated RMEs in Table fT|. 



The RMEs listed in Table |T| have been further multiplied by the factor 



327rag/ip 

where /i is the pd reduced mass and a is the fine structure constant, in order to match the 
different definitions used in the present study and by the authors of the fit. Here, X2J+1 
stands for either m2j+i, P2J+1 or q2j+i- Note that the phase of each RME is simply related 
to the elastic pd phase shift 6^^"^, as discussed in Refs. P,p2|. In particular, at these low 
energies to a very good approximation 6^^'^ ~ ctl, where ctl is the L-wave Coulomb phase 
shift. For example, the phase of the calculated doublet p2 RME at Ep = 35 keV is found 
to be 22.648°, which is to be compared with an elastic pd phase shift 5^^^ of 22.635" (at 
this energy, cri = 22.625°). As can be seen from Table [TITI , the most significant differences 
between theoretical and experimental RMEs are found for \p4\. 

The experimental value of the quartet to doublet ratio, tmi = 2|m4p/|m2p, for the Mi 
strength had been determined to be 0.49 ± 0.04 in Ref. 0] from an analysis of the T20 
data at 90°. The value obtained from the "experimental "Mi RMEs listed in Table |Tl| is 
0.43 ± 0.05. Both of these are in good agreement with the theoretical prediction of 0.475. 

It is interesting to analyze the ratio tei = |P4p/|p2p- Theory gives tei — 1, while from 
the fit it results that tei ~ 0.74 ± 0.04. As discussed above, the calculation of these RMEs 
is not infiuenced by uncertainties in the two-body currents, since their values are entirely 
given by the LWAc form of the Ei operator, which has no spin-dependence. It is therefore 
of interest to examine more closely the origin of the above discrepancy. If the interactions 
between the p and d clusters are switched off (hence reducing the pd scattering wave function 
^1+2 '^^^^ in Eq. ( p.4|) to the product of a deuteron wave function times Yi{rpd) Fi{qrpd)) , the 
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relation rEi — 1 then simply follows from angular momentum algebra, apart from negligible 
corrections due to the small P-wave components of the ^He wave function. Deviations of 
this ratio from one are therefore to be ascribed to differences induced by the interactions 
in the L,S,J = 1,1/2,1/2 and 1,1/2,3/2 wave functions. To study these differences, we 
define the "density functions" p2j+i{rpd) with the property 



P2J+1 = / drpdP2j+i{rpd) ■ (4.24) 

•/ 

In Fig. I, the functions ^[p2{rpd)] and ^[p4,{rpd)] are displayed with and without including in- 
tercluster interactions. As expected from the analysis above, the two functions ^[p2j+i{rpd)] 
are indistinguishable when these interactions are ignored. When the latter are included, the 
functions ^[p2j+iirpd)] are shown by the thick dashed (J = 1/2) and thick solid (J = 3/2) 
lines. Intercluster interactions have a significant effect reducing, however, both integrated 
values P2J+1 by the same amount ~ 10%, with the result that ^[p2] is still ~ ^[p^]- 

The AV18/UIX interactions produce essentially the same asymptotic behavior in the 
J =1/2 and J=3/2 doublet P-wave scattering states. This is directly confirmed by a com- 
parison of the calculated nuclear pd elastic phase shifts with the AV18/UIX Hamiltonian 
model in the and ^^3/2 channels, as shown for a few incident proton energies in Ta- 



ble The and ^^3/2 phase shifts are found to be very close to each other over the 

whole energy range considered in the table. Note that the AV18/UIX predictions for these 
phase shifts at Pp = 3 MeV are in excellent agreement with the values extracted from the 
phase-shift analysis (PSA) performed in Ref. 0. 

Therefore, the "experimental" value tei = 0.74 ± 0.04 is at variance with predictions 
based on the AV18/UIX Hamiltonian model. It should be emphasized that the present 
study ignores, in the continuum states, the effects arising from electromagnetic interactions 
beyond the static Coulomb interaction between protons. It is not clear whether the inclusion 
of these long-range interactions, in particular their spin-orbit component, could explain the 
observed splitting between the p2 and p4 RMEs. A calculation incorporating them in the 
continuum is currently underway, and it will be the subject of a forthcoming paper. 
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Finally, the calculated Ay and T20 analyzing powers at Ec.m. = 75 keV and 100 keV are 
compared with the data of Ref. [§] in Fig. |^. Note that the theoretical predictions have 
changed slightly with respect to the earlier results in Ref. 0], due to the improvement in 
the description of the pd P-waves discussed previously. Somewhat better (worse) agreement 
between theory and experiment is now found for T20 (Ay) than previously reported [Q. The 
large discrepancy at small angles for T20 is present also at these energies. 



B. The pd capture at 2 MeV 

Measurements of capture polarization observables at -Ec.m. = 2 MeV have been reported 
recently by Smith and Knutson who also extracted, by a fitting procedure, values of 
the contributing RMEs at this energy. To reduce the number of parameters in the fit, they 
made use of both pd elastic scattering and radiative capture data. In fact, as discussed in 
Ref. the phases of the radiative capture RMEs are fixed by the elastic S-matrix, when 



no channels other than elastic or radiative capture ones are open. Furthermore, by using 
the invariance of the nuclear Hamiltonian under parity and time-reversal transformations, 
it can be shown that the quantities |^ 

^LSJ ^ ^.,s^sj ^ uis,„s'Ee''' , = e-^'-'' E Uis,,'S'Mf''' , (4.25) 

L'S' L'S' 

are real. In the above expression 6^^"^ are the eigenphase shifts and Ulgj^,g, is the mixing 
matrix, into which the ^-matrix describing pd elastic scattering can be decomposed (both 
the eigenphases and mixing matrix elements are real under the deuteron breakup threshold). 
The RMEs {X = E, M) used in the expression above are related to the RMEs Xf^'^ 

defined in the present work (see the Appendix) and in Ref. @] by 



vLsj 

^LSJ ^ ^-if^s-j^p . (4.26) 

In Ref. the real quantities £ and M. have been used as free parameters and extracted 
from the data. The measured observables, though, were not sufficient to univocally select the 
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values of the RMEs and two sets of parameters have been presented 0. In order to compare 

directly with the quantities 8 and M., we have transformed our RMEs via Eqs. (|4.25|) 

and ( |4.26| ), using the eigenphase shifts and mixing matrices predicted by the AV18/UIX 

model. The RMEs obtained in this way are reported in Table [V| along with the two sets 

of "experimental "values. Inspection of the table indicates that significantly better overall 

agreement exists between "set 1 "and the theoretical RMEs. We will only consider this set 

in the discussion to follow. 

At this energy, the two doublet E'l's are the dominant RMEs. They have similar values, 

—ill — li^ 

and are in good agreement with those extracted from the fit. Note that E^^'^ E^'^'^ 
corresponds to \p2\ — \Pi\-, since the mixing induced by the matrix U in Eq. (|4.25| ) is negligible 
for these RMEs. Therefore, at this energy the relation |p2p/|p4p — 1 is well verified also by 
the "experimental "RMEs. 

The S-wave Mi RMEs (TW^^^ and A^^^^) are quite well reproduced by theory. Instead, 
significant differences between theory and experiment are found for the D-wave Mi RMEs- 
those extracted from the fit are an order of magnitude larger than predicted by theory. 
It should be emphasized, however, that most of the observables show sensitivity to these 
RMEs only in the small {9 < 30°) and large {9 > 150°) angle regions. However, using the 
experimental values for the D-wave Mi RMEs, rather than the calculated ones, does not 
produce "theoretical "observables in significantly better agreement with data. 

3 1 3 3 

The "experimental "quartet P-wave Ei RMEs {Si^^ and £^i^^) are overestimated by 
theory by almost a factor of 3. These RMEs are extremely sensitive to two-body currents. 



as can be seen in Table VI. The values obtained in the LWAc, LWAb and LWA2 approxima- 



tions are also listed in Table |V|. Similar considerations to those discussed in the previous 
subsection apply to these RMEs. Here we only point out that the LWAc values of the quar- 
tet El RMEs are rather close to corresponding values extracted from the data (set 1). This 
is particularly evident when the comparison is performed directly in the observables, as can 
be seen from Fig. ^. In the figure, the dotted and thin solid curves are obtained in lA and 
by using the full current, respectively, while the thick solid curves are obtained by retaining 
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the values calculated in the LWAc approximation for all electric dipole RMEs. 

Note that for the observables a, Ay, and iTn, the contributions of the quartet Eis are 

negligible. These observables depend mainly on the doublet Eis and S-wave Mi's, and 

are well predicted by theory. In contrast, the observables T20 and T21 depend linearly on 

the quartet -Ei's, and are not well reproduced by the "full" theory. The differences between 

the "full" and "experimental" values for these RMEs suggest that the present model for the 

two-body currents may have deficiencies. Note that the LWAc results for the electric dipole 
3 1 ^ 3 3 

transitions Si^'^ and Si^^ happen to be close to the "experimental" values, as shown in 



Table 0, and hence good agreement is obtained between the experimental and LWAc- 
calculated T20 and T21 observables. It is, however, important to emphasize that such an 
agreement is to be considered purely accidental, since the next to leading order contributions 
LWAb and LWA2 are comparable to the leading order LWAc results. 

Some of the differences found between the theoretical and experimental RMEs could be 



due to the different elastic eigenphase shift and mixing angle parameters, used in Eq. ( |4.25| ) 



to define the real RMEs S and Ai. To clarify this point, the elastic eigenphase shifts 
obtained with the AV18 /UIX model are compared with those used by the authors of Ref . 
in Table |V11|. There is good overall agreement between the two sets of parameters. We 



note that the small differences in the eigenphase shifts are responsible for the large 
underprediction of the Ay and iTn observables in elastic pd scattering. 

As discussed in the previous section, the LWAc and "full" estimates of the quartet Ei 
RMEs are rather different (see Table ||). However, when the LWAc values of these RMEs 
are used in evaluating the various observables, the resulting changes are found to be rather 
small, even for T20 and T21. In this energy regime, the observables are in fact dominated by 
the doublet Ei and the S-wave Mi RMEs, and the contributions of the quartet Ei RMEs are 
found to be negligible in all cases. At higher energies, though, the quartet Ei becomes of the 
same order of magnitude of the S-wave Mi RMEs, and therefore have large effects on some of 
the observables. Clearly, the significant differences between the "full" and LWAc predictions 
for the quartet Ei RMEs have a large impact, particularly for the tensor observables. 
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Finally, in Fig. ^ the "full" results for the S'-factor are compared to data in the cm. 
energy range 0-2 MeV stored at the web site [http : / /pntpm . ulb . ac . be/nacre . htir . 



C. Contributions to the Gerasimov-Drell-Hearn integral 

The Gerasimov-Drell-Hearn (GDH) sum rule connects the helicity structure of the photo- 
absorption cross section to the anomalous magnetic moment of the nuclear target, and is 
derived using Lorentz and gauge invariance, crossing symmetry, causality and unitarity of 
the forward Compton scattering amplitude p4| , |35| . For the case of ^He, it is given by 



^ r^jm^JLM = 2A. , (4,27) 



where (Tp^^(a;) are the photon absorption cross sections in which the photon helicity and 
^He spin are either parallel (P) or antiparallel (A), msHe and the ^He mass and 

anomalous magnetic moment, and c^th is the threshold energy. As discussed in the Appendix, 
it is related to the inclusive Rt' response, measured in polarized electron scattering from a 
polarized spin 1/2 target. 

In ^He, the GDH integral is h-ne = 498 /ib, using the experimental value Kane = —8.366 
for its anomalous magnetic moment. The photodisintegration threshold in ^He is u;th=5.495 
MeV, corresponding to pd breakup. It is useful to divide the integral into the part up to pion 
production threshold, and the part above this threshold. For the part above pion threshold, 
the ^He nucleus should have roughly the same strength as the neutron, i.e. ~ 230 yub. Such 
an expectation is based on the fact that the ^He ground state consists predominantly of a 
spherically symmetric S-wave component, in which the proton spin projections are opposite 
and the net polarization is therefore due entirely to the neutron. Ignoring corrections to this 
naive estimate, it is expected that the tu-region from the photodisintegration threshold up 
to the pion threshold should contribute about 266 /ib to the ^He GDH integral. A realistic 
description of the pd and ppn continuum states for energies above the three-body breakup 
threshold in terms of PHH wave functions, while certainly possible, is not yet presently 
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available. A calculation of the contribution to the GDH integral from the ppn threshold up 
to the pion threshold is therefore not possible. In the present work, however, we study this 
contribution in the 2 MeV window where only two-body photodisintegration channels are 
energetically allowed. 

In fact, the TUNL group has recently made the first experimental determination of the 
contribution to the GDH integral of ^He from the energy region up to 53 keV above the pd 
threshold |jl2|. In this region, there are only six dominant RMEs corresponding to electric 



and magnetic dipole transitions, as already noted before. Thus, ignoring the contributions 
from higher order multipoles, we find that the cross section difference in Eq. ( [4.27|) is simply 
given by 



IGn'^afip 



UJ 



"^4^ I |2 , I |2 , k4p 



-1^21 + - \P2\ + ^ - \Q2\ + 



(4.28) 



2 2 ' 2 

in the notation of Eqs. ([4.15|) -( ^.17|) . In particular, for the purpose of comparing with the 



discussion of subsection A, note the factor l/v2J+T difference between the X2j+i and X2j+i 
in Eq. (g^). 



As discussed in subsection A, the TUNL group has determined the relevant RMEs from 



an analysis of the polarized capture data, the details of which are reported in Ref. |]T2 
Because of time reversal invariance, the RMEs for the capture reaction are related to those 
for the photo-absorption reaction by phase factors, which are irrelevant for the Act defined 
above. 

It is convenient to define 

m = r du^jM^lM. , (4.29) 

with, obviously, I{uj oo) = hue- The experimental values obtained by integrating up 
to uJi = 5.522 MeV and uJ2 = 5.548 MeV are presented in Table |V1I1| , where they are 
compared to predictions obtained by including one-body only and both one- and two-body 
currents (columns labelled lA and FULL). Note that these values represent small negative 
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contributions to the total strength expected below pion threshold. Table [VII1| also lists the 
individual contributions to /(cJi) from the — |m2p+ |m4p/2, — |p2p+ b4p/2, —|q'2p+ 194^/2 
RME combinations (rows labeled Mi, Ei 5=1/2, and Ei 5=3/2, respectively). 

The total contributions including the Mi, Ei S=l/2 and Ei 5=3/2 strengths are found 
to be in lA an order of magnitude smaller (in absolute value) than data. This is because 
in lA |m2p — 0.5 x |m4p, |p2p — 0.5 x |p4p and the quartet (5=3/2) Ei strength is very 
small. 

The ratio ^0.5 for the doublet to quartet Mi strength obtained in lA is consistent 
with predictions for pd capture at zero relative energy obtained with the Faddeev method 
using a variety of realistic Hamiltonians [^] (the ratio is found to have only a weak en- 
ergy dependence). When two-body currents are included, the doublet Mi strength |m2p 
becomes roughly twice as large as the quartet Mi strength |m4p, a result also consistent 
with the earlier calculations [0. This makes the overall Mi contribution to I{uj) negative 
and relatively large (in absolute value). 

However, the nearly exact cancellation between the doublet Ei strengths |p2p and \p^\'^/2 
(or 1^2 P and |p4p) discussed earlier, is not significantly influenced by the inclusion of two- 
body currents. The quartet Ei RMEs remain negligible. Thus the total contribution to I{uJ) 
is mostly due to Mi strength. While the results in the "full" calculation are in much better 
agreement with data than those in lA, a factor of two discrepancy persists between theory 
and experiment. This discrepancy is mostly due to the difference between the calculated 



and measured 5=1/2 Ei strength, as the second row in Table |VIII| makes clear. It should 

also be noted that the measured I{uj2 = 5.548MeV) = — 1.120±0.218 nb is a very tiny piece 

of the expected contribution to the GDH integral of ^He below pion threshold, i.e. 266 /xb, 

and that it has the opposite sign. 

The result for /(tJa) with uj^=7.7 MeV corresponding to the threshold for complete 

breakup into ppn is 1.1 fih. We find that the cross section difference Aa changes sign in 

the "full" calculation at about uj=6 MeV. In the range of energy 6 MeV < uj < 7.7 MeV the 

1 1 J 

polarized cross section is dominated by the i?i^ RMEs. Indeed, the |p2p and |p4p strengths 
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are three orders of magnitude larger than the strength from any of the other contributing 
RMEs. However, the large cancellation between — |p2p and |p4p/2 persists also in this uj- 
region, although to a lesser extent than found above: 85% of Aa is due to these terms, while 
the remaining 15% comes from the other RMEs. Note that in this region, RMEs other than 
those included near the two-body breakup threshold need to be considered. For example, the 

2-- 

electric quadrupole ^ at = 7.5 MeV {Ep = 3 MeV) is found to give a 4% contribution 
to Aa. 

While the results reported above indicate that an extremely tiny piece of the total sum- 
rule strength is located in the threshold region, we find, by direct comparison with experi- 
ment, that this integral observable is very sensitive to the effects of two-body currents. The 
inclusion of these currents reduces the discrepancy between theory and experiment from a 
factor of ten to a factor of two. Further studies are needed to understand the physical origin 
of the difference in the leading P-wave Ei RMEs responsible for the remaining discrepancy. 



D. Threshold electrodisintegration of '^He 

The most recent and systematic experimental study of the threshold electrodisintegration 
of ^He and we are aware of was carried out by Retzlaff et al. [|l3l at the MIT/Bates Linear 



Accelerator Center. The longitudinal and transverse response functions and Rt were 
obtained using Rosenbluth separations for three-momentum transfers in the range 0.88- 
2.87 fm~^ and excitation energies from two-body thresholds up to 18 MeV. The ^H(e, e') 
data are the only measurements at these energy and momentum transfers. Inclusive ^He 
electron scattering data from earlier experiments [0 are in agreement with the Retzlaff et 



al. measurements, after scaling for the slightly different kinematics. 

The Rl and Rt ^He data at momentum transfer values g=0.88, 1.64 and 2.47 fm~^ are 
compared in Fig. ^ with calculations using PHH wave functions obtained for the AV18/UIX 
Hamiltonian model, and one-body only (dashed lines) or both one- and two-body (solid 
lines) charge and current operators. Note that the contributions associated with the L=0- 
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5 pd scattering states are retained in the calculation, which is then fully converged. No 
calculations of the response functions have been carried out at this time. There is 
satisfactory agreement between theory and experiment for all cases, but for the longitudinal 
response at g=2.47 fm~^. The two-body components of the electroexcitation operator play 
an important role, particularly for the transverse response at the highest g-values. The 
relative sign between the one- and two-body contributions is consistent with that expected 
from elastic form factor studies of ^He. There it is found that two-body current (charge) 
operators increase (decrease) the lA predictions for the magnetic (charge) form factor at 
g < 3 fm~^. 

A feature of the longitudinal ^He data is the presence of an enhancement near thresh- 
old. Such an enhancement is particularly pronounced at low q. It is not observed in the 
longitudinal response [jl3[. At the momentum transfers under consideration here, the 



longitudinal strength is almost entirely due to a Co transition involving the dominant S- 
wave component of ^He and the doublet S-wave pd scattering state. It has been shown by 
Heimback et al. that the enhancement results from the constructive interference of the 
amplitudes in the two-body breakup of ^He that correspond to the virtual photon coupling 
directly to a proton or to a correlated proton-neutron pair. In ^H, the virtual photon for 
the two-body breakup channel can couple only to the correlated pair, since coupling directly 
to a proton leaves an unbound neutron pair and thus a three-body final state. The present 
calculations correctly account for this threshold enhancement in the longitudinal response 
of ^He. 



Previous calculations, reported on in Ref. [0, used either the Faddeev equations ||39 



or the orthogonal-correlated-state (OCS) method to describe the bound and scattering 



wave functions. Both calculations did not include two-body charge and current operators. 
However, the former used the central Malfiiet-Tjon interaction, while the latter was based 



on a realistic Hamiltonian including the older Argonne f 14 two-nucleon [^] and Urbana 



VII three-nucleon interactions. The longitudinal and transverse response functions are 
surprisingly well predicted in the threshold region by the Faddeev calculation, but are both 
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underestimated, particularly the transverse response, in the OCS calculation. The latter 
also fails to reproduce the observed enhancement in longitudinal strength between the two- 
and three-body breakup thresholds, probably because of the approximate treatment of final- 
state-interaction effects between the proton and deuteron clusters. However, in view of the 
importance of two-body currents (see Fig. |^), the agreement between the data and Faddeev 
results is presumably accidental. 

Finally, in Figs. | and || we show the Rl, Rlt', Rt and Rt' response functions at a fixed 
excitation energy of 1 MeV above the pd threshold in the three-momentum transfer range 
0-5 fm~^. In Rl and Rlt' the S-wave pd continuum states give the dominant contribution, 
while in R^ and Rt' both S- and P-wave states give equally important contributions over 
the whole q range. All response functions are substantially affected by two-body currents, 
however, the sensitivity to these is particularly pronounced for Rlt' and Rt'- 

In Fig. |l^ we show the unpolarized cross section, and the Alt' and Ati asymmetries in 
the threshold region at an incident electron energy of 4 GeV. The asymmetries are relatively 
large at high g, and particularly sensitive to two-body currents. Note that Alt' has been 
found to be little influenced by uncertainties in the electric form factor of the neutron (the 
Galster parametrization is used in Figs. |7|^I]|), except at the highest g- values. The cross 
section for the chosen kinematics (incident electron energy of 4 GeV, fixed pd excitation 
energy of 1 MeV, and 0° < < 14°) is dominated by the longitudinal response function. 
Note that in Fig. [1^ we also show the plane-wave-impulse-approximation (PWIA) results. 
These have been calculated by approximating the wave function 

^>{li'^{VW\k)= [[s^®U^^)]s®YL{h)\J J''^^J''^ • (4-30) 

cyclic ijk ^ P^pd 

The lack of orthogonality between the ^f:^2"'"^' (PWIA) with L = S* = 1/2 and the ^He 
ground-state wave functions is responsible for the excess cross section at low q obtained in 
PWIA with respect to the "full" calculation. 
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V. CONCLUSIONS 



We have reported calculations of pd radiative capture observables at energies below the 
three-body breakup threshold, and of longitudinal, transverse and interference response 
functions measured in polarized electron scattering from polarized ^He in the threshold 
region for momentum transfers in the range 0-5 fm~^. These calculations have been based 



on the Argonne vis two-nucleon |rD| and Urbana-IX three- nucleon interactions, and 
have used accurate bound and continuum wave functions, obtained with the PHH method 
PJ3|JT^. The model for the electromagnetic operator has been taken to consist of one- 
and two-body components, the latter ones constructed consistently with the two-nucleon 
interaction [jTOipSf . In recent studies, this theory has been shown to correctly predict the 



static properties of the trinucleons pSf and A = Q nuclei as well as their associated 
elastic and transition electromagnetic form factors. 

A satisfactory description of all measured pd observables has emerged with the exception 
of the differential cross section and tensor analyzing power at small angles for Ep < 40 keV, 
and the contributions to the GDH integral at energies in the cm. range 0-53 keV. A com- 
parison between the calculated RMEs and those extracted from fits to the measured data 
has shown that the large p4 RME associated with the channel L, S, J = 1, 1/2, 3/2 is overes- 
timated by theory at very low energy. It has been speculated that at these energies (below 
< 50 keV) long-range electromagnetic interactions beyond the static Coulomb repulsion 
between protons might play a role. These electromagnetic interactions have already been 
included in bound-state calculations (such as those reported here), where they have been 
found to contribute to ground-state energy differences between mirror nuclei or members 
of isomultiplets IQ. However, they have yet to be included in scattering-state calculations. 
Work along these lines is currently underway. It should emphasized that at higher energies 
the splitting between the p2 and p4 RMEs appears to be much reduced, as the fits to the pd 
observables at 2 MeV indicate. This result is further corroborated by a phase shift analysis 
of pd elastic scattering at energies below the ppn breakup threshold [Q, which leads to very 
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close values for the and ^^3/2 phases. At 2 MeV the tensor observables T20 and T21 

are sensitive to the "small" quartet Ei RMEs. The "full" calculation fails to correctly predict 
these observables, suggesting that the present model for the electromagnetic transition op- 
erator may have deficiencies. Finally, the validity of the long-wavelength approximation has 
been analyzed. This approximation has been found to be inadequate for the calculation of 
inhibited (and hence "small") electric dipole transitions. Model calculations of low-energy 
capture and photodisintegration observables, based on the long-wavelength form of the Ei 
operator (the LWAc operator of Eq. (^78|) ), should therefore be viewed with suspicion unless 
explicitely verifying that higher order corrections are indeed negligible. 

The body of data available from ^He(e, e') inclusive scattering experiments in the thresh- 
old region is not as extensive, at present, as that from pd capture experiments. However, 
polarized electron scattering data from a polarized ^He target below the ppn threshold should 
become available in the near future |^^. The longitudinal and transverse data measured at 
Bates by the Saskatoon group are in reasonable agreement with theory, although the data 
have rather large errors. The crucial role played by two-body charge and current operators 
should be re-emphasized. 
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APPENDIX A: EXPRESSIONS OF OBSERVABLES IN TERMS OF REDUCED 

MATRIX ELEMENTS 



Expressions for the angular distributions of the cross section, vector and tensor analyzing 
powers and photon linear polarization coefficient were derived in terms of reduced matrix 
elements of electric and magnetic multipole operators in Ref. 0. It it useful to extend 
that analysis also to the case of the threshold photo- and electro-disintegration of ^He, 
since at the the small excitation energies of interest here only a relatively small number of 
electromagnetic multipoles are expected to contribute significantly. 

The electromagnetic transition amplitudes between an initial ^He bound state with spin 
projection (T3 and a final pd continuum state having proton and deuteron with relative 
momentum p and spin projections, respectively, (J2 and a, are given by: 

p<xa2<x3(p,q) = (^i,:iajp(q)l^3,i<x3) ' (Ai) 



p,q) = (vI/W |eA(q)-j(q)|vl/3i ) , (A2) 



where q is the momentum transfer and eA(q), A = ±1, are the transverse polarizations of 
the (real or virtual) photon. The wave function with ingoing-wave boundary condition is 
expanded as 

SSz LLzJJz 

where the '^^^2^''^ ^ are related to the '^1^2'^" introduced in Sec. || via 

L S 

Here is the Coulomb phase shift and is the i?-matrix. Introducing the expansion 
above into the matrix elements, Eqs. (|Al|)-( 1A2|) , one finds: 



jt,.3(P.q)=4vr E {-■^f{^crAa2\SS.,){SS^,LL-Mj;)YLLX^)3'i'LM) ' (^5) 



LLzS S z J J z 
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j'Sialq) = (<2''^^"^I^A(q) • j(q)|^3,i.3) > (A6) 



and similar expressions hold for the Pcro-2o-3(P; amplitudes. It is now convenient to take q 
as defining the z-axis, i.e. the spin-quantization axis. Standard techniques |^5| then lead to 
the following expansions in terms of reduced matrix elements of Coulomb (C), electric (E) 

^f/3(q) andji'f^l 



and magnetic (M) multipoles for the amplitudes Pj^JJq) and jj^xcr,,{fi) 



^ °o Of \ ] 1 

J^ZM) = -^Eiy^7^(2^3,a| JJ.) [XM^''{q) + Ef^'^iq)] . (A8) 

The matrix elements Pj^J.^ (qz) and jj^x„,_^ (qz) are calculated with the Monte Carlo inte- 
gration techniques discussed in Ref. 0, and from these the reduced matrix elements C^^-^i^q), 
Mj-^^q) and E^^\q) are obtained via Eqs. ([A^)- (|A8|) , for example 

Ct'hq) = ^pfkq^) , (A9) 

VZTT 2 2 



M^'^q) 



2 2 '->'-'-> 



(AlO) 



'2t^ 

The inclusive cross section for polarized electron scattering from a polarized spin 1/2 
target can simply be written as |^ 



j:{q,u) + hA{q,u) , (All) 



d^duj 

uj) = aMlvLRdq, tu) + vrRriq, ^^)] , (A12) 



A{q,u!) = (7 M[vLT'RLT'iq,^^) sin 9* cos (p* + VT'Rx'iq,^) cos 9*] , (A13) 

where aM is the Mott cross section, the coefficients Va are functions of the electron kinematic 
variables, h = ±1 is the helicity of the incident electron, and the angles 9* and (p* specify 
the direction of the target polarization with respect to q. The response functions Ra con- 
tain the nuclear structure information. We note that the sum over the three-nucleon final 
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states, implicit in their definition, is restricted to include only the pd continuum, since the 
excitation energies of interest here are below the threshold for the three-body breakup. The 
longitudinal-transverse and transverse-transverse asymmetries Alt' and At' are defined as: 

VLT'RLT'{q,Uj) 



ALT'{q:<^) 
AT'{q,L0) 



VLRL{q,u;) +VTRT{q,uj) ' 

VT'Rx'iq, ^) 



(A14) 



VLRL{q,uj) +VTRT{q,uj) 

Explicit expressions for the response functions Ra in terms of reduced matrix elements 
of electromagnetic multipole operators are easily obtained: 



Rl — fpd XI 

LSM 



Rt = fpd {\E^^'^\^ 
LSje 



LSJ\2 



(A15) 



(A16) 



LSJ 



2J+1 



( c!i'' + i ^ ) * [/7^ ( m!:'' + e!i'' ) 



-i ( M^^^ + E^^^ ) 



, (A17) 



-2 J(J+3/2)(J-l/2)a 



(A18) 



where the phase-space factor fpd is given by fpd = and in the interference response func- 
tions the notation X^^"^ for the reduced matrix elements means X^Jj^^^^- The magnitude 
of the relative momentum p is fixed by energy conservation 



Lo + — E2 -\- 



g2 2 

+ 7r ^ 



(A19) 



2(m2 -l- m) 2jj, 

where E2 and E^ are the two- and three-body ground-state energies, m2 is the deuteron 
mass and is the 1+2 reduced mass. 

The photo-disintegration cross section is simply related to the Rt response function: 
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a/(uj) = ^^Rt{uj) , (A20) 

UJ 

where for real photons q = uj. In particular, the difference of cross sections in the integrand 
of Eg . ([4.271) is easily related to the response function Rt' 

aVcu) - aViu) = ^^Rt,{uj) . (A21) 

CO 
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TABLES 







a2 (fm) 




04 (fm) 




Th 


Exp 


Th 


Exp 


nd 


0.63 


0.65 ± 0.04 


6.33 


6.35 ± 0.02 


pd 


-0.02 




13.7 





TABLE L Predictions obtained from the AV18/UIX Hamiltonian model with the PHH method 
for the nd and pd doublet and quartet scattering lengths 02 and 04. 



lA FULL LWAc LWAb LWA2 

e\1^ -24.2 -26.2 -26.0 -0.1 0.1 

e\P -6.1 0.6 2.1 0.6 -0.4 

e\1^ 33.1 36.7 36.9 0.0 0.1 

e\1^ -2.7 0.1 0.6 -1.3 0.7 

TABLE II. Doublet and quartet Ei RMEs in fm^/^ calculated with the AV18/UIX Hamil- 
tonian model for the reaction ^H(p,7)^IIe at zero energy in lA and FULL approximations. The 
contributions LWAc, LWAb, and LWA2 are reported in columns 4-6, see text for an explanation. 
Note that the LWA2 contribution has been calculated by using only the magnetization term of the 
single-nucleon current, i.e. by only retaining the second term on the r.h.s. of Eq. ( |3.6| ). Statistical 
errors associated with the Monte Carlo integrations are in the range 1-5%. 
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RME 


lA 


FULL 


FIT 




0.172 


0.322 


0.340±0.010 


7714 


0.174 


0.157 


0.157±0.007 


\P2\ 


0.346 


0.371 


0.363±0.014 




0.343 


0.378 


0.312±0.009 


TABLE III. 


Magnitudes of the leading 


: Ml and Ei RMEs for pd capture 


at Ep = 40 keV. The 


values listed in 


the fourth column have b 


een determined in Ref. [|l^ from 


a fit to the measured 


observables, while those listed in the second and third columns have been obtained in calculations 


using either one 


-body (lA) only or both one- and two-body (FULL) currents. Note that the RMEs 


listed in this table are adimensional. 






^p(MeV) 






^-^3/2 


0.035 




-0.00392 


-0.00391 


0.213 




-0.699 


-0.695 


2.0 




-4.89 


-4.84 


3.0 




-7.37 


-7.17 



3.0 (PSA) -7.41(0.08) -7.18(0.04) 



TABLE IV. Nuclear elastic pd phase shifts (in degrees) obtained for a few selected proton 
energies with the AV18/UIX Hamiltonian model. The values extracted from the phase-shift analysis 
(PSA) of Ref. m are listed in the last row. 
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xliVillj 


set i 


set z 


AV io/ UlA 


nil 


n 001 
— U.zzi 




-O.z4o 


93 1 
y v I -1 




-0.355 


0.023 


0^^ 


n 001 
— U.zzi 


— U.zzO 


n 1 no 
-O.lyz 


9I 3 
y 1 1 


0.272 


0.476 


0.018 


933 
A /( 2 2 


U.io4 


U.oz4 


u.uiy 


-lii 
c 2 2 

^1 


z. / zi 


AO A 


z.o ^8 


_l3 1 


-0.122 


-0.118 


-0.308 


_ll3 
S['' 


2.742 


2.717 


2.771 


1 3 3 

£['' 


0.080 


0.061 


0.175 


q3 3 


0.061 


0.085 


0.072 



TABLE V. Real RMEs (xlO^) for pd radiative capture at E^.m. = 2 MeV, see text for def- 



initions. Sets labelled 1 and 2 are the RMEs obtained in Ref. from a fit to the measured 
observables. The RMEs in the column labelled AV18/UIX have been obtained by the calculation 
presented in this paper, including the full current operator. Note that the RMEs defined in this 
table are adimensional. 
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RME 


set 1 


lA 


FULL 


LWAc 


LWAb 


LWA2 




2.721 


2.352 


2.678 


2.699 


0.012 


-0.021 


1 3 1 

£\^^ 


-0.122 


-1.166 


-0.308 


-0.127 


0.071 


-0.050 


_ll 3 


2.742 


2.395 


2.771 


2.716 


-0.042 


0.011 


1 3 3 


0.080 


0.405 


0.175 


0.111 


0.121 


-0.069 


q3 3 


0.061 


0.069 


0.072 


0.070 


0.002 


0.001 


TABLE VI. Real RMEs 


(xlO^) for pd radiative capture 


at Eq ycl 


= 2 MeV as in 


Table 0. Set 



1 are the RMEs obtained in Ref. from a fit to the measured observables. The RMEs in the 
columns labelled lA, FULL, LWAc, LWAb, and LWA2 have been obtained in calculations using 
various approximations for the Ei operator as in Table 0. Note that the RMEs defined in this 
table are adimensional. 
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VV C*i V V' 


5'^'^"' Exp 


:LSJ AV18/UIX 




V o/ 




2q 


-24.9 


-27.9 


^1/2 


4.27 


4.27 


^ Sq /o 


116.1 


116.9 


-'-^3/2 


9.83 


9.98 


^^3/2 


4.12 


4.01 


^Pl/2 
1/2 


-1.82 


-2.15 


-'^1/2 


26.9 


27.4 


^^3/2 


-1.95 


-1.95 


'^^3/2 


29.4 


29.3 


"^-^3/2 


10.4 


10.4 



TABLE VII. Eigenphase shifts 6^^'^ for pd elastic scattering at Ec.m. = 2 MeV. The values of 
LSJ are given in column 1 in the format "^^^^Lj. The experimental phase shifts are those from 
Ref. Q, whereas the theoretical ones are calculated with the AV18/UIX interaction model. The 
values reported here are the phase shifts induced in the wave functions by the nuclear plus Coulomb 
potential. 
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cu (MeV) 


m 


FIT 


lA 


FULL 


5.522 


Ml 


-0.0524 ± 0.0077 


-0.0029 


-0.0609 


5.522 


El S =1/2 


-0.0361 ± 0.0095 


-0.0030 


+0.0027 


5.522 


El S =3/2 


-0.0026 ± 0.0007 


-0.0050 


-0.0001 


5.522 


Total 


-0.0911 ±0.0123 


-0.0109 


-0.0583 


5.548 


Total 


-1.120 ±0.218 


-0.161 


-0.582 



TABLE VIII. The contributions /(cJ) (in nb) for two energies uJ. In the third column, the 
values obtained from a fit of the experimental pd capture data are reported ||T^. The results of 
the theoretical calculations using the AV18/UIX Hamiltonian model and either one-body only or 
both one- and two-body currents are listed in the fourth and fifth columns, labelled lA and FULL 
respectively. The lines denoted by Mi, Ei 5=1/2 and Ei 5=3/2 report the partial contributions 
to /(w=5.522 MeV) of the corresponding RMEs. 
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FIGURES 

FIG. 1. The 5-factor for the ^II(p,7)^He reaction, obtained with the AV18/UIX Hamiltonian 
model and one-body only (dashed line) or both one- and two-body (solid line) currents, is compared 
with the experimental values of Refs. [0j8|. 

FIG. 2. The energy integrated cross section a{9)/ao {Anao is the total cross section), vector 
analyzing power Ay {6), tensor analyzing power T2o{9) and photon linear polarization coefficient 
P-yiO) obtained with the AV18/UIX Hamiltonian model and one-body only (dashed line) or both 
one- and two-body (solid line) currents are compared with the experimental results of Ref. 0]. 

FIG. 3. The functions ^[p2j+i{rpd)] calculated with the AV18/UIX Hamiltonian model for 
the states J = 1/2 and 3/2. The functions obtained by switching off the nuclear pd intercluster 
interactions are displayed by the thin dashed and solid lines for the J = 1/2 and J = 3/2 scattering 
states, respectively. The two lines are indistinguishable. When the nuclear interactions between 
the d and p clusters are taken into account, the functions are shown by the thick dashed and solid 
lines for the J = 1/2 and J = 3/2 scattering states, respectively. 

FIG. 4. Proton vector analyzing power Ay and deuteron tensor analyzing power T20 for pd 
capture at Ec.m. = 75 and 100 keV, obtained with the AV18/UIX Hamiltonian model and one-body 
only (dashed lines) or both one- and two-body currents (solid lines). The experimental values are 
from Ref. [^. 

FIG. 5. Differential cross section, proton vector analyzing power, and the four deuteron tensor 
analyzing powers for pd capture at Ec.m. = 2 MeV, obtained with the AV18/UIX Hamiltonian 
model and one-body only (dashed lines) or both one- and two-body currents (thin solid lines), 
are compared with the experimental values of Ref. Q. The results obtained in the approximation 
LWAc for the Ei operator are also shown (solid lines). 
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FIG. 6. The S'-factor for the ^H(p,7)^He reaction, in the cm. energy range — 2 MeV, obtained 
with the AV18/UIX Hamiltonian model and one- and two-body currents (soUd hne) is compared 
with the experimental values listed in the web site http : //pntpm . ulb . ac . be/nacre . htni| . 



FIG. 7. The longitudinal and transverse response functions of ^He, obtained with the 
AV18/UIX Hamiltonian model and one-body only (dashed lines) or both one- and two-body (solid 



lines) charge and current operators, are compared with the data of Ref. |13] at excitation energies 
below the ppn breakup threshold. 

FIG. 8. The longitudinal {Rl) and longitudinal-transverse {Rlt') response functions of '^He, 
obtained with the AV18/UIX Hamiltonian model and one-body only (thick dashed lines) or both 
one- and two-body (thick solid lines) charge and current operators, are displayed at a fixed exci- 
tation energy of 1 MeV for three- momentum transfers in the range 0-5 fm^^. The contributions 
associated with the (dominant) S-wave pd scattering states are also shown. 

FIG. 9. Same as in Fig. ^, but for the transverse (Rt) and transverse-transverse {Rt') response 
functions of ^He. The contributions associated with both S- and P-wave pd scattering states are 
also shown. 

FIG. 10. The inclusive cross section, and the vIlt' and Ax' asymmetries, obtained with the 
AV18/UIX Hamiltonian model and one-body only (dashed lines) or both one- and two-body (solid 
lines) charge and current operators, are displayed for ^He at a fixed excitation energy of 1 MeV 
for three-momentum transfers in the range 0-5 fm^^. The results in PWIA (dotted lines) are also 
shown. The incident electron energy is 4 GeV, and the electron scattering angle is in the range 
0-14°. 
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